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Abstract - We introduce an analytically solvable model for a fragmented object 
that, despite of a low degree of randomness and of the extreme simplicity of the 
breaking process, displays non self-averaging effects in its thermodynamic limit. 



In this paper we present an exactly solvable model for a fragmented object 
that, despite of the extreme simplicity of the breaking process, displays non 
self-averaging effects in its thermodynamic limit. The presence of non self- 
averaging quantities is a striking feature of a variety of statistical mechanical 
models of disordered systems, originating in condensed matter physics, as the 
low temperature spin glasses [Q, as well as population biology 0, dynamical 
systems' theory and mathematics |Q. In most cases a complete analytic 
study is not possible, and the interesting probability densities are obtained 
numerically. The model we propose is sensibly less random than all of the 
above models, and yet, quite surprisingly, has the same statistical properties. 

The question of self-averaging arises in disordered models when one con- 
siders the statistical fluctuations of a given thermodynamic or statistical ex- 
tensive property X of the model. Suppose to have a certain system of size 
(under "size" we shall understand a linear dimension of the system as well as 
the number of parts forming it, or equally a measure of its phase space, such 
as the number of its possible configurations) in which disorder is represented 
by some quenched random variables. For finite A^, X is sample-dependent, 
in the sense that to each sample of the system, namely to each realization of 
the quenched disorder, corresponds a unique value of X, and the ensemble 
average {X) of X is obtained by averaging over all possible realizations of the 
quenched disorder. The sample-to-sample fluctuations of X are described by 
the quantity 
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which clearly depends on and is non-zero for finite A^. If ^^{X) ^ in 
the thermodynamic limit N —>■ oo then X is said to be self-averaging and 
a sufficiently large sample is a good representative of the whole ensemble 
(think of the energy of a classical ideal gas of particles). But if V^{X) 
tends to a finite positive value, then X remains sample-dependent even in 
the thermodynamic limit, and an evaluation of X made on a no matter how 
large sample is not significant of the value of X on other samples. If this is 
the case, X is said to be non self-averaging. 

What one finds in many disordered models, some of which we cited above, 
is that in the thermodynamic limit the system's phase space appears broken 
into infinite basins, and that: (i) such breaking is sample dependent; (ii) for 
each breaking there are finite sized basins; (iii) the sizes of the basins are non 
self-averaging quantities. To see this, one usually considers the probability 
Y that two randomly chosen configurations (points in phase space) belong 
to the same basin. 



where Wg is the probability that a randomly chosen configuration belongs to 
the s-th basin, and the sum is made over all basins. To prove (ii) one has to 
show that y's ensemble average, calculated over all possible ways in which 
the phase space may be broken, is finite and positive, as is explained in Ref. 
IQ. If one then manages to show that Y lacks of self- averaging, (iii) would 
soon follow. Finally, (i) is a mere consequence of (iii). Hence one tries to 
demonstrate that 

• (Y) > (for feature (ii)); 

• var(F) > (for features (iii) and (i)). 

The most interesting quantity to compute is the probabihty density n(F) of 
Y, such that Il{Y)dY gives the probability that the value of Y for a random 
sample lies between Y and Y + dY. The fact that Y has a finite variance 
in the thermodynamic limit tells us that Il{Y) remains "broad", and this is 
confirmed by computer simulations, but an analytic derivation of 11(1") is 
not possible in most models. Along with n(y), other interesting quantities 
are the probability densities of the weights of the largest, second largest, . . . 
pieces in each sample, that are expected to remain "broad" as well, and the 
average number of pieces of a given size for a sample. 

The "broken phase space" feature explains why fragmentation models 
may attract some interest. One assumes that an object of size 1 is broken 
into infinite pieces by means of a certain process, which depends on a number 




(2) 



s 



2 



of quenched random variables. A realization of the breaking process is given 
by a specific sampling of the random variables and the result of each breaking 
is an infinite set of pieces of sizes Wg {s = 1,2,...). Averages are then 
calculated over all possible ways of breaking the object. Of course, not all 
breaking processes lead to non self-averaging effects. Suppose for example 
that the object is broken uniformly into M basins of size Wg ~ each; then 
y ~ jj, and both Wg and Y would go to zero in the M ^ oo limit. In some 
sense, the fact that a uniformly broken object is self- averaging suggests that 
a non-uniformly broken object, such as the randomly broken object P|, may 
be regarded as the simplest complex system. We shall now see that a much 
simpler, though similar, breaking process than that of the randomly broken 
object leads to comparable non self- averaging effects, with the advantage of 
being analytically solvable. 

We consider an object of size 1 and suppose to break it into infinite pieces 
by means of the following process: given a real number p G [0, 1], we tear the 
object in two pieces of sizes 1 — p and p respectively; then we do the same 
thing with the piece of size p, obtaining two pieces of sizes (1 — p)p and p'^ 
respectively, plus the one of size 1 — p obtained at the first breaking step. 
If we repeat the same procedure with the pieces of sizes p"^, p^, ... that are 
obtained at the second, third, . . . breaking steps, keeping the pieces of sizes 
(1 — p)p, (1 — p)p^, . . . obtained at the same steps we finally have a set of 
pieces of sizes 

Wi = 1 - p 
W2 = {l-p)p 

(3) 

Wg = 
. . . , 

where Wg denotes the size of the piece kept at the s-th step. Clearly, Wg = 
— p)p'^~^ = 1- We call such process a geometrical breaking, since the 
sizes of the resulting pieces form a geometric sequence. Now suppose that 
p is a random variable with probability density p{p). In this case we can 
imagine that a breaking sample is produced by choosing a random value of 
p from p{p), so that averaging over all samples means simply averaging over 
all possible values of p. 

As usual, let Y = J2g value for a single sample is given by 

y = 5:(i-p)V(-^) = il^ = l^, (4) 
^ 1 — p^ 1 + p 
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and the ensemble average of Y is simply 

(y) = \^p{v)dp. (5) 

When p is uniform we easily obtain the value 

(F) = log4- 1 ~ 0.386... . (6) 

The probability density n(y) over the samples may be calculated from 
the relation 'n.{Y)dY = p{p)dp expressing the conservation of probability. 
From formula @ we get p = {1 + Y)~^{1 — Y) and thus 



n(F) 

namely 

n(y) = 



dp 



dY 



Y 



[1 + Y)^'^\1 + Y 



y's ensemble average is given by (Y) = YIl{Y)dY and one can verify that 
for a uniform p the value (|^) is recovered. It suffices to use the fact that 

2Y ( \ \ 

-(iF = 2 log|l + y| + - - + constant. (9) 



The second moment (F^) = /J y2n(r)dr may also be calculated in a 
straightforward manner for a uniform p. We can make use of the relation 

to obtain (F^) = 3 - log 16 ^ .227. . . 7^ (y)^. The variance of Y is finally 
given by 

var(F) = (y^) - ^ 0.078 ... . (11) 

The fact that Y has a non zero variance in the thermodynamic limit, that 
for this model is represented by the infinite pieces in which the object is 
broken, proves that F, and consequently the sizes VFg, lack of self-averaging 
for a geometrically broken object. We remind that all numerical values were 
obtained under the assumption of a uniform density p{j)). 
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Let us now turn to the other interesting quantities. Since for each sample 
one may find a piece of maximum size, then a second maximum sized piece 
and so on, one might be interested in the densities of the probabihty that the 
n-th largest piece in a sample has size between W and W+dW. Let us denote 
such densities by Pn{W), so that Pi{W)dW represents the probability that 
the largest piece in a sample has size between W and W + dW, P2{W)dW 
represents the probability that the second largest piece in a sample has size 
between W and W + dW and so on. For a geometrically broken object one 
can easily understand that, for each sample, 

Wi>W2> ■■■ (12) 

since, for p e [0, 1], it is always 1 — p > (1 — p)p > (1 — p)p^ > ■ ■ ■ . Hence, 
the largest size is 1^ = 1 — p, so that, from the relation Pi{W)dW = pip) dp 
and from the fact that, for the largest piece, we have p— 1 — W, follows 



Pi{W) 



dp 



dW 



p{l-W)^ p{l-W). 



(13) 



For the second largest piece one has W = {1 — p)p- According to a standard 
method of probability theory we can write 



P2{W) = 



p(Ph 



+ 



p{p. 



\W'{p+)\ \W'{p 

where p+ and p_ are the roots of the second order polynomial (1 —p)p — W 

dW 
dp 

1 



(14) 



0, and W'{p) = ^ = 1 - 2p. We have 



which leads to 



P2{W) 



1 



p±^-(l±^l^^), 



(15) 



(16) 



In principle, one can calculate the densities PniW) for every n by solving 
the polynomial (1 — p)p^^^^ — W = and applying the technique used above 
in its generalization given by the formula 



^|M^'(p«)|' 



(17) 
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where p^^^ is the i-th root of the polynomial. 

The average largest, second largest, . . . sizes are obtained by {Wn) 
J WPn{W)dW which may be as well written 



For a uniform p this leads to 



f ^p{,p)dp. 
Jo 



{Wn) = --^ = -r^, (19) 
n n + 1 n\n + Ij 

(n = 1 corresponds to the average largest size, n = 2 to the average second 
largest size, and so on), so that on the average the largest piece is three 
times as big as the second largest, six times bigger than the third largest, 
. . . , fortyfive times as big as the 9-th largest, and so on[|. 

We can also calculate the average number of pieces of size between W 
and W + dW , denoted by f{W). We know that the n-th largest piece has 
size Wn = (1 — p)p"'~^. This size is maximum for p = jo^™'™) = Hul, in which 
case it is given by 

^(max) _ ^^^^(max)) _ I ) _ ^20) 

This means that each Pn{W) is defined just for W < ivi""^"^ (in fact, P2{W) 
is defined for W < j = 14^2™^^^). If in a certain sample we have a piece of 
size ^ < W < 1 it must be the largest piece in the sample, since the second 
largest piece can have a maximum size of |. Hence, for ^ < < 1 one has 

f(W) = P^{W) = p{l-W), (21) 



by means of (13). If we have a piece of size ^ < W < j it can be either 



the largest one or the second largest one, since the third largest one has a 
maximum size of ^. Hence, for ^ < < | one has 

f{W) = Pi{W) + P2{W) 



p{l~W) ^ 



VI -41^ 

H.p(i(i_yr3w)) 



f,(^^(i + vr^)) + (22) 



^ This type of distribution of the average weights is known as Zipf's law (M. Gell-Mann, 
Zipf 's law and related mysteries, lecture held at the SFI Complex Systems Winter School, 
Tucson, January 12-24, 1992 (unpublished)). 
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by means of (|13|) and (|T^). In the same way, if we have a piece of size 
(n+i)"'+^ < < — (for any positive integer n) this can be any of the 
first n largest pieces, since the {n + l)-th largest piece can have a maximum 
size of , ■ Hence for , ,^^+1 < W < — one has 

f(W) = Pi{W) + --- + Pn{W). (23) 

This shows how one can reconstruct the function f{W) in all intervals of the 
form <W < , and hence for all W e]0, 1]. 

To summarize, we have shown that this geometrically broken object has 
the same statistical properties as other, well known disordered models, and 
displays lack of self-averaging. Furthermore, the simplicity of its definition 
makes it possible to solve it exactly, and we have derived all the most inter- 
esting quantities, such as the probability densities n(F), Pi{W) and P2{W), 
and the average number of pieces of size W, fiW). We believe that this 
model is interesting because the breaking process involves just one random 
variable, whereas in the randomly broken object the number of random vari- 
ables entering the breaking process is infinite. In this sense, the geometrically 
broken object is less disordered than other disordered models, and the fact 
that comparable non self-averaging effects are present is quite surprising. 

We have applied the results obtained in this paper to the study of a 
stochastic evolutionary model P], in which the phase space breaks like a 
geometrically broken object, to show for example that the model lacks of 
self-averaging in the thermodynamic limit. 
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